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SUMMARY 


This series of reports relates the status of research on the fabrication of 
glancing incidence optical components at the Goddard Space Flight Center. This 
work has been performed by the Optics Branch of the Experimental Fabrication 
and Engineering Division under SRT Program 283-188-38-01-10. 

The investigation has four aspects — theory, design, fabrication and testing. 
This report is primarily concerned with the elementary theory and design of a 
two element glancing incidence telescope. The GSFC X-ray Telescope Ex- 
periment for ATM, a Type 1 telescope (concave parabola followed by a concave 
hyperbola), is analyzed using a modified ray tracing program outlined by 
Wolter \ It is shown that the image quality of such a system can be significantly 
improved by matching the detector to the optimum focal surface and by vignet- 
ting rays which contribute to the flare in comatic images. Vignetting is accom- 
plished by reducing the hyperbola element to approximately one half of its 
length from an unvignetted configuration. The total weight of the hyperbola ele- 
ment is, therefore, reduced by one half, resulting in a savings of material cost 
and fabrication time. The imaging properties of a one -third scale ATM type 
telescope are discussed. 


iii 


X-RAY OPTICS 


PART I 


INTRODUCTION 

The desire of scientist to extend observations of primarily solar phenomena 
into the spectral region below 500 A has spurred the investigation of using 
glancing incidence type optics to obtain spatial as well as spectral information 2,3 
The gathering of spatial information in the X-ray region is not new as some in- 
vestigators have used an array of devices ranging from pinhole cameras to vari- 
ous mirror arrangements 4 ’ 7 . However, surfaces of rotation appear to have 
unique characteristics worth investigation. These are the component rigidity 
and relatively large photon gathering capacity combined with the ability to pro- 
vide spatial information. 

An excellent paper by Wolter 1,4 describes in detail some of the fundamental 
properties of systems which are under investigation at GSFC. Wolter, however, 
approached the problem from the aspect of applying these optical systems to X- 
ray microscopy. Furthermore, it has been evident that a more detailed account 
of the imaging properties is necessary to acquaint the scientist with the theoreti- 
cal limitations of such systems and also to reveal to the optical technician what 
one may expect to see when fabricating non -normal incidence components and 
systems. 

It is important to note that the term "non-normal incidence" is taken to be 
synonomous with "glancing incidence" for it is precisely from this non-normal 
incidence characteristic that imaging problems arise. 

The principal reason for using non-normal incidence configuration is to 
realize the efficient reflection of energy in the X-ray region. Figure 2 illus- 
trates this fact quite clearly, for as the angle of incidence proceeds from normal 
incidence to near normal incidence and finally to non-normal incidence the 
amount of energy reflected increases. 

A parabola of revolution might appear as the first likely candidate for use 
as a glancing incidence element. However, objects imaged off axis by this con- 
tour suffer severe aberrations. It has been demonstrated that the reason for 
this effect is that imaging characteristics of single glancing incidence elements 
grossly violate the Abbe Sine Condition *. This may be understood if one consi- 
ders that a perfect image of a point source is formed at the center of curvature 
of a collapsing spherical wavefront. Therefore, it is the function of the optical 
component to convert the wavefront emanating from the object point into a 
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spherical wavefront which will collapse onto the corresponding image point. The 
Sine Condition is written in the form 8 

h 

— : 7T ~~ constant - r (1) 

sin a ' ' 

for a point source at infinity, where h is the radial height from the optical axis 
at which the incident ray strikes the aperture and 9 is the angle the reflected 
ray makes with the optical axis. Any deviation from the principal spherical 
surface defined by Equation (1) may then be thought of as contributing to imper- 
fections in the image. Clearly, a single glancing incidence element does not 
conform to the Sine Condition. 

However, for a combination of glancing elements such as the concave 
parabola-concave hyperbola (labeled type 1 by Wolter) illustrated in Figure 2, 
the noncompliance of the principal surface with Equation (1) or more precisely 
the offense against the Sine Condition is given bv 1 


OSC = Ar 


n 



( 2 ) 


The above combination may then be designed to minimize OSC and also serve as 
an efficient X-ray collector. If the effective aperture is given by A, then 


where 


S “ r A a 2 

Furthermore = 40 where the maximum slope angle 6 is small 

2. m a x xr o max 

(i.e., 6 ~1°) then 

' 7 max ' 


A r a 

y 

J p m l n 

Hence, to minimize OSC for a given minimal collecting area A, and maximum 
glancing angle 6 max, (equal to the maximum slope angle of the parabola for an 
on axis ray) the diameter of the aperture must be increased. 
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Designs for type 1 telescopes are rather inflexible as all of the definitive 
parameters are fixed if the maximum slope angle, minimal collection area and 
focal length f = r , are specified. From Appendix I and Figure 2 


yp min 

f sin 46 - y . . 

max ^nmin 

(3) 

P 

= ^min tan 0 max 

(4) 

C 

= f/2 

(5) 

a 

= f (2 cos 26 - 1 )/2 

v max 7 

(6) 

b 

- ( c 2 “ a 2 )' 4 

(7) 


TA , T* 


y 

J p m a x 

- + y 2 . 

77 J p m l n 

(8) 


The radial image size, r , in the Gaussian focal plane can be calculated for 
a half field angle 8 from the equation 


T 


f sin S 


1 - cos 40 

max 

2 cos 46 

max 


However, the above formula only applies with accuracy for imaging by the 
zone of intersection of the parabola and hyperbola. It is therefore necessary to 
develop a ray tracing procedure for a more comprehensive assessment of the 
imaging properties of the telescope. 


RAY TRACING PROCEDURE 

A set of 19 rays are traced from points on a circle defined by the intersec- 
tion of the YZ plane with the parabolic aperture. Rays are traced from points 
on this circle spaced at ten degree circumferential increments of an angle 4> 
about the optical axis, starting with the positive y axis (<£ = 0°) through 4> - 180° 
to the negative y axis. It is necessary to only trace rays through one half the 
aperture because of the symmetry of the telescope. A total of eleven sets 
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(denoted as SI, S2, ... Sll,) or 209 rays are traced for each half field angle, 5 . 
The X intercept of each set of rays is such that ten percent of the full aperture 
area is encompassed in the zone defined by two successive sets of rays. The 
locations of the YZ ray intercept or focal plane is moved in and out of the 
Gaussian focus (M = 0). Positive values of M indicate that the focal plane is 
moved toward the telescope whereas negative values of M indicate that the focal 
plane is moved away from the telescope. Half field angles were varied from 
S = .000 to S = .006 radians in .001 radian increments. 

Details of the ray trace procedure which closely follow that outlined by 
Wolter, are outlined in the Appendix II. A 7094 computer was used to perform 
the computations. 


TELESCOPE DESIGN 

The first telescope design to be traced using the program described above 
is that of a type 1 (concave parabola followed by a concave hyperbola) configura- 
tion. Using the formulation of Equations 3, 4, 5, 6, 7 the configuration of the 
telescope is given by the following parameters if f = 75 inch, 6 max = 0.916 
degrees and A = 2.3 in 2 

p = 0.07663156 inch y . = 4.792896 inch 

a = 37.4616644 inch y pinax = 4.868791 inch 

b = 1.6951988 inch 
c = 37.500000 inch 

This design was used as the basic configuration for the GSFC ATM X-ray 
telescope experiment. 


IMAGING PROPERTIES 

The spot diagrams in Figure 3A depict the image of a star formed by the 
telescope in the Gaussian focal plane for various half field angles. The image 
becomes severely elongated in the XY or tangential plane as greater field 
angles are viewed. As the focal plane is moved toward the telescope, M > 0, 
the image size for large values of S decrease in size, however, for large dis- 
placements of the focal plane (M > .020) there is corresponding increase in the 
images made at small field angles, Figure 3B. 
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The magnitude of this effect on the tangential spot diameter versus focal 
plane location is shown in Figure 4. As is evident in the diagram, there is a 
leveling off or a best focus region for each half field angle. This arises from 
the fact that the telescope suffers from curvature of field, i.e., the tangential 
rays incident on the front and rear aperture of the parabola are not imaged in 
the same plane. The longitudinal and laterial displacement of these rays results 
in an image which appears to be sheared as illustrated by the half aperture spot 
diagram in Figure 5. It is evident that although the image can be decreased in 
size by displacement of the focal plane toward the telescope that the downward 
or comatic flare of the image is still present to some degree. 

One may use a vignetting technique to minimize the flare effect, first by 
locating where the tangential rays are incident on the hyperbolic surface, 

Figure 6, correlating this information with their relative displacement in the 
XY plane, and then reducing the length of the hyperbola such that the unwanted 
rays are vignetted. The procedure for using this technique is as follows: The 
hyperbola X intercepts are plotted versus circumferential angle 4> for each half 
j field angle S and every set of rays. There is but one X intercept on the hyper- 
j bola for each of the 11 sets of rays for S = 0. For half field angles greater 
I than 8 - 0, the X intercepts vary along the length of the hyperbola within a 
J given set and have common cross over positions approximately located at 4> = 
90° and 4> = 270°. It is seen that for field angles greater than 8=0, using the 
8 = 0, X intercept as a reference that the X intercept within a given set of 

rays is closer to the front of the telescope for circumferential angles where 
90° < 4> £ 270° and correspondingly the X intercept moves away from the front 
of the telescope and thereby elongates the hyperbola for 0° £ < 90° and 

270° < 4> < 360°. 

As stated previously, the tangential rays from the upper and lower rim of 
the parabola do not have a common focus. It is seen from Figure 7 that the 
meridonal tangential rays, y p max , focus is closer to the telescope than that of 
the lower tangential rim rays, y p min - There are therefore two cross over loca- 
tions of these rays between the two tangential foci. (It is in between these cross 
over points where the leveling off or minimization in image size occurs as was 
noted in Figure 4.) The ray contributing to the downward most portion of the 
image comatic flare in the YZ plane of the crossover point closest to the 
Gaussian focal plane arises from the front zone of the parabola, y p m . ix , at 
,6 = 0°. The cross over point is approached as £ - 180°. It is therefore evi- 
dent that those rays incident on the meridonal zone on the parabola in the 
vicinity of <f> - 0° are the source of the comatic flare. However, it is precisely 
these rays which are also incident on the last portion of the hyperbola. Hence, 
if a back portion of the hyperbolic element is eliminated, these rays will be 
vignetted and the comatic flare suppressed to some degree depending on just 


5 



how short the hyperbola is made. In order to minimize the vertical diameter of 
the image, the length of the hyperbola was selected such that the last ray incident 
on the lower rim of the hyperbola for 4 > = 0°, intersects the ray originating from 
4 > - 180° at y p mi n in the same YZ plane which contains the cross over point 
closest to the Gaussian focal plane. This is shown in Figure 7 for S = .004 
radian. For the design under investigation, this technique results in halving the 
length of the hyperbola. 

The loss of energy" in the image (assuming a one-to-one correspondence to 
the decrease in effective aperture) is plotted in Figure 8 versus the length of 
the hyperbola. Vignetting at large field angles is substantial, but a considerable 
reduction of image size is realized. A comparison of unvignetted and vignetted 
images, Figure 9, for each field angle illustrates the image compacting quite 
vividly. The angular subtense of the tangential diameter of these images in 
their optimum focal positions is indicated in Figure 10. 

The shape of the optimum focal plane can be approximated as a hyperbola 
described by the equation 


(x + a ) 2 y ' 2 

a 2 ” b 2 


1 


where 

a = .08476 inch 

b = .30151 inch 

However, it is not always possible to shape the image detector to the optimum 
focal surface. It is therefore informative to examine the distribution of energy 
within an image in a flat focal surface as well as the optimal focal shape. As 
stated previously, eleven zones of the parabola were traced, each pair of zones 
encompassing 10% of the field aperture area. The distribution of energy within 
an image may then be correlated to the area circumscribed by these sets of 
rays. This energy (area) distribution was measured by incrementally summing 
the distribution of the area for each zone in the ± Z direction within the YZ focal 
plane. There is an obvious increase in energy" at the image core, Figure 11, as 
one progresses from the Gaussian to optimal focus. (The arrows on the Y axis 
indicate the point at which vignetting takes place if the hyperbola is reduced to 
one half of its length from the unvignetted configuration.) 
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Resolving Power 


To what degree this telescope will resolve two closely space point sources 
is difficult to predict. It is clear that by using selective vignetting and shaping 
the detector to the contour of the optimum focal surface that the relative percent 
of the light in the hard core of the image is increased. This judgment has been 
made strictly on the basis of geometrical ray tracing. However, there are other 
factors which must be considered when predicting image quality. 

First, one must realize that in practice, the resolving power will also be 
dependent upon the optical quality of the contour used for imaging. Geometrical 
optics can validly be applied to a surface where the wavelength is small com- 
pared to the dimensions of the surface roughness. In this case, the reflectance 
in a given direction would only be a function of the slope of the surface imper- 
fection. As the wavelength becomes longer diffraction and scattering become 
important factors and the reflectance is determined by the size and spatial fre- 
quency as well as the slope of the imperfections. Finally, as the wavelength be- 
comes large compared to the dimensions of the surface roughness, diffraction 
phenomena dominates the reflectance characteristics of the surface. Unfortu- 
nately, the above considerations are further complicated when perfect conductors 
are not used as the reflecting element and one must therefore be concerned 
about polarization as well as phase and amplitude effects. 

To further illustrate the magnitude of the problem for predicting the image 
quality of an x-ray telescope, consider the data in Table I as applied to the tele- 
scope designed to form images in the spectral region 6A to 26A. From Table I 
it is clear that the roughness of optically polished surfaces span the spectral 
region of interest. Variations in the measurements of the root mean square 
roughness, a, of a given material are attributed to melt differences and the 
techniques applied to polish and measure the surface roughness. Because it is 
desired to make a surface as smooth as possible it is preferable to use and 
polish a material such as fused silica or Cer-Vit using the bowl feed technique. 
Semi-conductors, conductors, and conductive coatings have not been considered 
here because of the lack of quantitative data. But, even for the super smooth 
surfaces of fused silica or Cer-Vit, the wavelength is not large compared to the 
surface roughness. It is therefore to be expected that diffraction and scattering 
effects will be convolved with the geometrical imaging characteristics of soft 
x-ray telescopes. 

The maximum permissible surface roughness height may be calculated from 
a strictly geometric viewpoint as illustrated in Figure 12. The optical path dif- 
ference, S, between rays A and B is calculated to be given by: 



Table I 


Roughness of Optically Polished Surfaces 


Material 

a 9 ~ 14 

r m s 

Flint glass 

154A (a) , 41A (f) 

Plate glass 

68A< a \25A (d >, 40A/ e \ 18A (f) 

Glass Slide 

56A< b > 

Pyrex 

51A< a \ 14A (f) 

Fused silica 

36A (a) , 3A (d) , 9.6A (f) , 2.7A (d \ 7A (d) < 

Borosilicate crown glass 

10A (c) 

Cer-Vit 

6.6A (d >, 6.9A/ d \ 10.lA< d \ 7A< d > < 


(a) Fresh feed polishing technique 

(b) Fire polished 

(c) Not stated (Ref. 11) 

Bowl feed polishing technique 

(e) Not stated (Ref. 10) 

(f) Not stated (Ref. 9) 


S 


h 


/ 1 - cos 2 9 \ 
\ sin 6 ) 


where 6 is the glancing angle. One finds by applying the Rayleigh limit to the 
optical path difference that 


— ( sin §. \ 

4 \ 1 - cos 2 6 ) 


(9) 
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Therefore, at a glancing angle of 6 = 1°, h /X ~ 7.1 which, for X. - 25A, 
corresponds to approximately one-thirtieth of a wavelength in the visible. This 
analysis, however, treats only the phase relationship of the focused beam. 

The effects of surface roughness on the amplitude of the reflected beam are 
quite complicated. The most thorough research in the literature 9, 15 has been 
for cases of normal incidence and X/o- << 1. 

A mathematical treatment of the subject by Porteus 16 has given some in- 
sight as to the complexity of the problem. The reflected beam is considered to 
be divided into two parts, namely, that which is coherently reflected and that 
which is incoherently reflected. The light which is incoherently reflected gives 
rise to wave fronts which are uncorrelated in phase unless they originate from 
a localized region on the surface, whereas the coherently reflected light consists 
of mutually interfering wave fronts from all parts of the reflecting surface. The 
coherently reflected component, R c , is related to the limiting reflectance R Q 
(i.e., as cr - 0) by the expression 



/• + CO 


2 


COS 2 y 


D ( z) e [ -47Ti (z/A) cos V] dz 


(10) 


where is the angle of incidence, z is the height of the surface irregularity 
and D(z) is the marginal density function of the surfaces heights which is related 
to the total surface joint density function Dj (s, z , z' ) by 


D(z) 



z' ) dz 1 


Although the surface heights z and z' are statistical variables at the points 7 
and 7', respectively, they are not mutually independent because of the continuity 
of the surface and therefore maintain a correlation which depends on 7 - 7' 
represented as s. 

The incoherent term for the reflectance, R it from a uniformly rough iso- 
tropic surface is given by 
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( 11 ) 


R 

1 

~r~ " x 




z' ) 


D(z)D(z' ) j e ^ ^ (Z } 1 d S ^dzdz' dsdf) 


where K d represents the propagation vector of the diffracted wave fronts. 

Porteus has shown that although Equations (10) and (11) may be used to describe 
the reflectance characteristics of rough surfaces at normal incidence where 
X > > a, there still exists a great deal of uncertainty when trying to predict 
short wavelength behavior of the incoherent reflectance (i.e., as X - cr). The 
reason for this difficulty arises from the fact that it is demonstrated that obvious 
parameters such as the root mean square roughness cr, the surface slope m and 
the autocovariance length a do not necessarily govern incoherence reflectance 
in this range because of their common derivation from the autocovariance 
function 


< Z Z* > “ <z> < z' > 


It is further demonstrated that subtle differences in the choice of a joint 
density function markedly influence the short wave length behavior. The reason 
that the above complexities have been emphasized is as follows. 

Consider an isotropic surface consisting of facets of random size and shape 
which are aligned parallel to a reference plane but at random levels relative to 
this plane. Further choose that the probability, p(s), that two points chosen at 
random on the surface a distance apart will lie at an equal level relative to the 
reference plane. If then p(s) and D(z) are required to be Gaussian the coherent 
reflectance for normal incidence is given as 


R c " R o eX P [” (477aA) 2 ] 


while the incoherent reflectance is 


R 


R o -^1 - exp [- ( 4 7T a/k ) 2 ] ^ -^1 - exp [- (-naa/X) 2 ]^- 
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where a is the semi vertex angle of the cone of acceptance. These last two 
equations describe very well the behavior of light reflected at normal incidence 
where a A < < 1. 

If a cos is substituted for a to account for the angle of incidence then 

R c = R o exp [ - (4770- cos \j//\) 2 ] (12) 


From the previous consideration of applying the Rayleigh limit it was shown 
for = 90° - 6 = 89° at 6 = 1° that cr/\ ~ 7.1. However, applying these values 
to Equation (12) it is found that 


Clearly a surface which closely approximates the joint density function 
described above and a roughness on the order of 7 \ is unsatisfactory. If the 
surface roughness is on the order of magnitude of the wavelength, o- ~ \ , such 

o 

as could be the case at k = 6A for a bowl feed polished fused silica surface 
then, 


R 


C 


R 


O 


'■v 

''V 


95 % 


This reflectance is far more acceptable than the previous example; however, 
it implies that the parabola and hyperbola contours can be polished to an accu- 
racy of A./500 in the visible region. It is more probable that one would achieve 
a surface which has a super smooth surface with a contour which only approxi- 
mates the aspheric contour. Then one has two surface roughness parameters, 
a l and a 2 , which describe respectively the super smooth surface and the 
"waviness" or larger variations from the prescribed contour. This latter con- 
dition is given as 


= exp £- (47 tA) 2 (o- 2 + O- 2 ) cos 2 


and 
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if one again assumes that cr can be replaced by a cos v . It is apparent that in 
order to achieve a reasonable coherent reflectance that a 2 - a l . However, 
again it is emphasized that the above treatment and experimental verifications 
have been for the case of normal incidence. There is a definite need for a pre- 
cise study of glancing incidence conditions. 

Secondly, the annular aperture characteristic of x-ray optics in theory 
should contribute to diffraction phenomena which would increase the theoretical 
resolving power of a diffraction limited system by narrowing the half intensity 
width of the principal maxima. However, the preceding ray trace results indi- 
cate that the system would only be diffraction limited over a field encompassing 
a few seconds of arc. Indeed, the on-axis performance is the only probable area 
where one would expect to see diffraction limited performance. A drawback of 
the annular aperture is that the redistribution of energy is such that the relative 
intensity of the secondary diffraction rings is increased when compared to the 
principal maximum as shown in Figure 13. The result is that although the Ray- 
leigh criterion can be exceeded, the background is also increased and a loss in 
contrast is to be expected. 

In summary, it should be noted that incoherent scattering, diffraction, and 
the geometrical imaging properties of the Type I soft x-ray telescope all tend 
to decrease the image contrast. 

Experimental Results 


Some interesting effects have been observed in the imaging properties of 
an x-ray telescope fabricated in this laboratory. The telescope, which was 
fabricated from borosilicate crown glass, is approximately one-third the scale 
of the GSFC ATM system described above. The ratio of the radii which define 
the annular aperture is defined as 
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y 

J p m i n 

6 s — = 0.984 

y pmax 

where 2 y p max = 82.6262 mm. The focal length, f , of the telescope was meas- 
ured to be 627.761 mm. 

The photographs in Figures 14A and 14B illustrate the telescope image 
quality when viewing a 0.58 arc second source in green light (\ = .5461^). If 
the concentric rings for the on-axis image arise from a diffraction limited 
quality, then the radial intensity distribution in the image must comply with the 
equation describing the diffraction by an annular aperture. This equation is 17 


I (P) _ 


(1 - e 2 ) : 


2Ji (x) J, (e x) 

e 2 2 


1 2 


e x 


(13) 


where 

2tt 

x “ ~ y p m a x 0 

<t> = angular ring radius at a point, P 
I (p) = intensity at point, P 
I o = incident intensity 

Jj(x) & Jj(ex) = Bessel functions of the first order 

The intensity minima are determined by setting the above equation to zero and 
solving the resulting equation 


Ji (x) - £ J, (ex) 


(14) 


Table II is a listing of x for minimas when e = 0.984. 

If the radii of the diffraction minimas as measured off the photographs are 
designated as r i , where i is the minima number counting from the center of 
the diffraction pattern, and M is denoted as the magnification incurred in the 
photographic process, then 


r £ = f 0. M 
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Table II 

Listing of x For Minimas When e = 0.984. 


Minima # 

X 

1 

2.436 

2 

5.566 

3 

8.722 

4 

11.889 

5 

15.052 

6 

18.216 


However, 




x. k 


2 tt y 


p m a x 


Therefore, 


2 TT 


p m a x 


X 


f X M 


r . 


k r 


(15) 


The quantity x. was calculated from Equation (15) using the system parame- 
ters outlined above and the measurements of when M = 43. The results are 
listed in Table III for k = 176.089 cm - 1 . 

It is clear from Table HI that the on-axis performance of the telescope is 
diffraction limited. The discrepancies in the larger values of x £ are in part 
due to the slightly oval shape of the ring pattern, which indicates that there is a 
defect in the pinhole used as a source and/or in the circular symmetry of the 
telescope optics. 
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Table III 


Minima # 

Minima Radii 

Experimental 

Theoretical 

Error 

i 

r . 

i 

X. 

l 

X 

% 

i 

.013000 cm 

2.289 

2.436 

-6.0 

2 

.033225 cm 

5.851 

5.566 

+5.1 

3 

.051350 cm 

9.042 

8.722 

+3.7 

4 

.071100 cm 

12.520 

11.889 

+5.3 

5 

.090375 cm 

15.914 

15.052 

+5.7 

6 

.108140 cm 

19.042 

18.216 

+4.5 


At 0.5° off axis the diffraction rings are still present; however, some flare 
in the image is evident as well as scattering. Figures 15A and 15B illustrate 
the loss in sharpness of the image on axis as the angular subtense of the object 
is increased to 2.32 and 5.79 arc seconds, respectively. 
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Figure 1. Reflection 




18 


Figure 2. Type 1 X-Ray Telescope Geometry 
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Figure 3A. Spot Diagrams at Gaussian Focus, M - 0.00 inch 
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Figure 3B. Spot Diagrams Inside Gaussian Focus, M - +0.05 inch 
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Figure 4. Tangential Spot Diameter Versus Focal Plane Locatic 
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Figure 5. Half Aperture Spot Diagram 
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Figure 6A. Hyperbolic Intercept Versus Aperture Circumferential Angle 
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Figure 6 C. Hyperbolic Intercept Versus Aperture Circumferential Angle 
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Figure 6D. Hyperbolic Intercept Versus Aperture Circumferential Angle 
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Figure 6F. Hyperbolic Intercept Versus Aperture Circumferential Angle 
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Figure 9A. Spot Diagram in Optimal Focal Position for Vignetted and Unvignetted Telescope 
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Figure 9B. Spot Diagram in Optimal Focal Position for Vignetted and Unvignetted Telescope 
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Figure 9C. Spot Diagram in Optimal Focal Position for Vignetted and Unvignetted Telescope 
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Figure 9D. Spot Diagram in Optima! Focal Position for Vignetted and Unvignetted Telescope 
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Figure 9E. Spot Diagram in Optimal Focal Position for Vignetted and Unvignetted Telescope 
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Figure 9F. Spot Diagram in Optimal Focal Position for Vignetted and Unvignetted Telescope 
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Figure 10. Optimal Focal Plane Locations of Vignetted Telescopes 
for Various Half Field Angles 
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APPENDIX I 


MATHEMATICAL DERIFICATION OF THE 
FUNDAMENTAL PARAMETERS OF A 
TYPE I X-RAY TELESCOPE 


The basic parameters of a Type I x-ray telescope (concave parabola fol- 
lowed by a concave hyperbola) are derived using the notation in Figure 1. Rays 
A and B parallel to the optical axis x are incident at maximum and minimum 
points, respectively, defining the extent of the parabola aperture. The dashed 
lines which intersect the points <x p , y p ) and [(x p min ,. y p min ) or (x h max , y h max )] 
are drawn as tangent lines to the parabola and hyperbola with slope angles 6 
and 6 ax , respectively. 

If then one is given the maximum glancing angle 0 max for an incoming ray, 
the focal length f and the require collecting area A , the following parameters 
may then be calculated: 

a, b, C, p, y p min , Yh max’ X hmax’ Ypmax’ X pmax’ '"min’ min’ X hmin 

where the parabola is given by 

y p 2 = p ( 2 % + p) 

and the hyperbola by 


( x h 



b 2 


1 


For y : v f sin 4 0, 

^pmin •'pmin m 


For p : P ' y pmin tan 


^p min 

For x . : tan 2 <9 = 

x pmin max X 

p rain 


x - 2 f 2 cos 2 2 6 

pmin m 
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For 

For 

For 

For 

For 

For 

For 

For 


a - 


r p m i n 


sin 2 8 sin 4 6 

max max 


a - — ( 2 cos 2 9 “1) 

9 v max ' 


- |"x - f cos 4 8 1 

2l_pmin maxj 


C 


b : b - (c 2 -a 2 )^ 


X ! X — X 

h max ‘ hmax p min 


max ' max ^p min 


y : A 

J p m a x 


77 


fy 2 - y 2 ■ ) 

V^pmax ^ p m 1 n / 


p m a x 


= r*- y . . i' 

|_ 77 J P m i n J 


X : V 

pmax ^ p m a x 


p( 2 x + p) 

^ v pmax ^ J 


pmax 


*t 


pmax 


“ P 


8 . : tan 0 

min min 


pmax 
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For y hm 


For x hm 


tan 


- 1 


p m a x 


y h 


p m a x 


C ± 


p m a x 


p m a x 


r p m a x 


(1 - b 2 ) + a 2 


b 2 


p m a x 


p m a x 


h min 


p m a x 


min 


p m a x 
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APPENDIX II 


ELEMENTS OF THE RAY TRACE PROCEDURE FOR TYPE I 
AND TYPE II X-RAY TELESCOPES 


Rays incident on the circumference of the parabolic element are considered 
to lie in the XY plane and make an angle delta, S, with the XZ plane as shown in 
Figure 1. 

The location of the point of intersection P (x, y, z) of the ray with the 
aperture is given by the angle phi, <t>, as measured from the Y axis. Before 
tracing the ray through the system, the point P (x, y, z) is rotated through the 
angle until it is coincident with the point P (x'j, y x , Zj) on the Y axis. The 
angles formed by the incident ray with the XY plane and XZ plane are denoted 
as S s and § t , respectively. These angles, § s and S , are calculated as shown 
in Figures 4 and 5. 

It is necessary to determine the above angles since the direction cosines of 
the ray striking p ( x j, y j, Zj ) are calculated as follows from Figures 2 and 3. 

For the case of 0° < <p < 90° 


COS 

a - 

COS 

II 

/— N 

CO 

o 

o 

00 

»— < 

“ cos 8 

cos 

yS = 

cos 

(90° - s t ) = 

sin S t 

cos 

7 = 

cos 

(90° + S S ) - 

- sin 8 


and for the case of 90° < . <fi < 180° 

cos a = cos (180° + S) = - cos S 

cos /3 = cos (90° + § t ) = - sin 8 t 

cos y - cos (90° + § s ) = - sin § s 


Figures 4 and 5 have been drawn for the cases when the point p (x, y, z) is 
located in the quadrant 0° £ <t> £ 90° and 90° < <fi < 180°. It is seen from Fig- 
ure 4 that for 0° _< 4> < 90° 
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sin 8 


AB 


N sin S 


N 


sin 8 


AC 


N 


sin cp 


AC 


AB 


o°o sin 8^ - sin S sin 4> 


sin S , 


AD 


N 


cos p - 


AD 


AB 


sin S t - sin 8 cos <p 


and from Figure 5 $ hat for 90° < <p < 180° 

AB 


sin 8 


N 


sin 8 


AD 


N 


sin (180° - <£) 


AD 


AB 


So sin 8 s - sin 8 sin 4> 


sin 8, 


AC 


N 


n , BD _ AC 

cos (180 <p ) Ab Ar 


<g, sin 8 - -sin 8 cos <p 


AC - N sin S s 

AC = AB sin <p - N sin 8 

AD = N sin S t 

AD = AB cos 4> - N sin 8 

AB = N sin 8 
AD = N sin S s 

AD = AB sin <$> - N< sin 8 

AC = N sin S t 

AC = “AB cos <p = - N sin 


sin <p 


cos <p 


sin <p 


8 cos <p 


It should be noted that sin S t may be considered in the form sin *> t - 
sin 3 cos cp since cos $ is positive for 0° < d> <90° and negative for 90° < < 

180°. 


The next parameters to be determined are the direction cosines of the 
surface normal to the point P (x’j, y x , z x ). The surface normal is denoted by 
N x as shown in Figure 6. The slope of the tangent to the point P (x' lf y 1? 
is defined as a /2. Calculation of the direction cosines is then as follows: 


cos a 



sin 



x 



cos /3 



- cos 



y 

Ni 


Hence: 


x 


sin — 


y = -Nj cos — 

The direction cosines of the surface normal are then: 


sin ~ 2 ~ . ~ cos — , 


0 


if N is considered to be a unit vector. It is clear from Figure 6 that a x may 
be calculated from the relationship 


a 


l 



where 
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y i - [p( 2 *'i + p)] a 

Direction cosines defining the reflected vector, Sj,a.re then calculated by 
using the equation: 


Si = S o - 2Nj (S 0 N 1 ) 

where the direction cosines have been shown to be S (-cos 8, sin 8 , -sin 8 ) 

/ a i a i \° 

for the incident vector and N t I sin — , - cos 0 J for the surface normal. 
Substitution of these directional components into Sj yields 


-cos 8 i + sin 8. j - sin 8 k - 2 

t J S 


- 


"1 _ 


sin — i - cos — j + O k 


x 


(- cos 8 i + sin § t j - sin S g k) 



cos 


a i _ _\ 

~2 j + 0 k j 


After vector multiplication and collecting terms 
Sj = (sin § t sin a x - cos 8 cos a x ) i 

+ (“ sin 8 t cos a l - cos S sin a x ) j - sin S s k 

The direction cosine’s defining the reflected vector are then denoted as 

S lx = sin S t sin a l - cos S cos a x 
S ly = -sin S t cos a x - cos 8 sin a x 
S- - - sin 8 

1 z S 

This reflected ray is then incident at point P ( x' 2> y 2 , z 2 ) on the hyperbolic 
surface which has the form, in the present coordinate system, of 
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( x 2 " c > 2 


1 


a 


2 


b 2 


b 2 


For the case of the Type I telescope, a concave parabolic mirror followed 
by a concave hyperbolic mirror, the above equation is correct. However, for 
the Type H telescope, a concave parabolic mirror followed by a convex hyper- 
bolic mirror, the hyperbolic surface is described by 


( x 2 + c ) 2 y 2 


z 

17 


1 


It is convenient to perform the linear transformation 


x 


l 


c 


X 


2 


X, 


c 


since now the hyperbolic element of both types of telescopes is described by 



1 


The computer program must consider the value c as multiplied by -1 in 
order to perform the proper transformation for the Type II telescope. 


Now it is necessary to calculate, 1 2 , the magnitude of the vector between 
points P (x l9 y l9 z x ) and P (x 2 , y 2 , z 2 ) as illustrated in Figure 7. The angle 9 
is used as an arbitrary supplementary angle in deriving the following equations: 


S lx - cos a - cos (180° + 9) - - cos 9 



x 


2 


*. + t, s u 
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cos yS 


cos (90° + 6) 


- sin 8 


S 


ly 




2 


) 


y 2 


Vi 


+ £ 


s 


ly 


S lz - cos y - cos (90° + 8) - - sin 8 



z 


2 


s 


1 z 


Substitution of these equalities into the hyperbolic surface equation has the 
form 


(*1 + ^2 Six) 2 _ (Vl + ' 2 S ly ) 2 _ «2 S 1 Z > : 

a 2 b 2 b 2 


By expanding and collecting terms one has 



\2 1 


/ b \ 2 

[(7 

s , 2 , - 1 

- 2 T 

y i s iy _ It) x i s ix_ 


+ [(t) x ' " y2 * b2 ] = 0 

or 

B 2 ± (B 2 2 - A 2 Q 2 )' a - root for Type I 

£ 2 = 1 

a 2 + root for Type II 


where 
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1 


A 


2 



B 


2 


Vl S ly 



s 


lx 


c 2 




- b 2 


The next step is to determine the direction cosines of the normal to the 
hyperbolic surface. The quantity, N 2 , is defined as the magnitude of the surface 
normal between point P (x 2 , y 2 , z 2 ) and where it intersects the X axis. This 
axial intersection takes place because the hyperbola is a surface of revolution. 
The X intercept is calculated from Figure 8 as follows. The equation for the 
hyperbola is written in the form 


1^( X 2> y 2 * z 2 ) 


2 -> 2 
x 2 y 2 z 2 

X = 0 

a 2 b 2 b 2 


since the equation of the normal may be determined by the general relationship 


x - x 2 y - y 2 z - z 2 

3F ( x 2 , y 2 , z 2 ) — ^ F( x 2 , y 2 , z 2 ) ~~ ^F(x 2 , y 2 * z 2 ^ 

3x 2 dy 2 d z 2 


if 


F(x 2 , y 2 , z 2 ) - 0 


Therefore, 
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JLt = ZZ. 

3x 2 g 2 

BF _ 2 Vj. 

^y 2 -b 2 

BF _ 2z 2 
Bz 2 - b 2 

and the equation of the normal to P (x 2 , y 2 , z 2 ) is determined to be 



The X intercept of the surface normal with the optical axis is thus found to be 


x 



y - 0 


The magnitude of N 2 is given by 


N 2 
1N 2 



z 


2 

2 


however 


y 


2 

2 



- b 2 


and by substitution 
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n 2 



T 


1 


N 


2 


The direction cosines of N 2 are then determined from Figure 8 and once 
again 6 is used as an arbitrary supplemental angle. 


N 2x - cos a - cos (90° - 6 ) 


sin 6 



N 


2 x 



T 


N 2y - cos /3 - cos (180° - 6) - cos 6 


y 2 


N 


2 


N 2y = - y 2 T 


N 2z - cos 7 


0 - z, 

N„ 


N 


2 z 


T 


The direction cosines of the ray reflected from the hyperbola are again 
determined from the relationship 

S 2 = S, - 2N 2 (S, N 2 ) 
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which becomes 


S lx 1 + S ly 


i + Si, k - 2 (N 2x i + N 2 j + N 2z k) 


[(Si, 


i + S ly j + S,. k) (N, v i + N 2v j t N 


2 x 


2y 


2 z 


E)] 


S 2 = S ly j + S„ k - 2 (N 2x i + N 2y j + N„ k) DEL 

where 

DEL = s i,N 2i * S ly N 2y t S„N 2 , 

S 2 = <S lK - 2N lx DEL) T + (S ly - 2N 2y DEL) T + (S„ - 2N 2 , DEL) k 

The direction cosines of the ray reflected from P (x 2 , y 2 , z 2 ) are now 
determined as 

S 2, = S u -2N 2 „DEL 

S 2, * S ly -2N 2y DEL 

s 2, = S lt - 2N 2 ,DEL 


Finally it is necessary to determine where this reflected ray intersects a 
plane perpendicular to the optical axis of the telescope. The linear transfor- 
mation which was applied above has moved the location of the Gaussian focal 
plane to the X intercept x = c. Hence, the X coordinate of the intercept or 
"picture" plane is 


x p - c + m 


where m is the amount this plane is displaced in the X direction from the 
Gaussian focus. 
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The magnitude of the vector connecting points P (x 2 , y 2 , z 2 ) and P (x F , y F , 
z p ) is denoted as T 2 and is illustrated in Figure 9. It then follows that 


S 2x - cos a - cos ( 180° + 6) - cos 9 



Hence 


X F ” X 2 + T 2 S 2x " c + m 




m 


-) 


and 


S 2y - cos /3 - cos (90° + 9) cos 6 

y F _ y 2 + ^2 ^2y 

S 2z = cos y = cos (90° + 6) - cos 6 - 

Z F ~ Z 2 + ^2 ^2z 

The true "picture" plane intercept P (x^., y‘ F , z' F ) is finally determined by 
reversing the original angular transformation through . This is accomplished 
by the familiar equations 

y F = y F cos 4> - z F sin 4> 

z'p = y F sin 4> + z F cos <£ 
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and 


Xp - c + m 
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